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Abstract 

The quantum Szilard engine (QSZE) is a conceptual quantum engine for understanding the 
fundamental physics of quantum thermodynamics and information physics. We generalize QSZE to 
arbitrary spin case, i.e. spin QSZE (SQSZE) and systematically study the basic physical properties 
of both fermion and boson SQSZEs in low temperature approximation. We give the analytic 
formulation of the total work. For the fermion SQSZE, the work might be absorbed from the 
environment, and the change rate of the work with temperature exhibits periodicity and even-odd 
oscillation, which is a generalization of spinless QSZE. It is interesting that the average absorbed 
work oscillates regularly and periodically in large-number limit, which implies that the average 
absorbed work in fermion SQSZE is neither an intensive quantity nor an extensive quantity. The 
phase diagrams of both fermion and boson SQSZEs gives the SQSZE doing positive or negative 
work in the parameter space of the temperature and the particle number of system, but they 
have different behaviors because the spin degrees of fermion and boson play different roles in 
their configuration states and corresponding statistical properties. The critical temperature of 
phase transition depends sensitively on the particle number. By using the Landauer’s erasure 
principle, we give the erasure work in a thermodynamic cycle and define a new efficiency (we 
refer it as information-work efficiency) to measure the engine’s ability of utilizing information to 
extract work. We also give the conditions under which the maximum extracted work and highest 
information-work efficiencies for fermion and boson SQSZEs can be achieved. 
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FIG. 1. The working cycle of the classical multi-particle Szilard engine, (a) A wall is isothermally 
inserted into the box. (b) Then a measurement is performed to know the number of particles of 
both sides, (c) The system does work by isothermally moving the wall to the equilibrium position, 
(d) The wall is removed isothermally and the system is restored to the initial state. 


I. INTRODUCTION 

The Szilard engine (SZE) is a single molecule engine designed for solving the Maxwell’s 
demon paradox[? ]. As is shown in Fig. 1, N ideal identical particles are prepared in a 
one-dimensional box of length L, and the system is always in contact with a heat bath of 
temperature T. Then an adiabatic wall is isothermally inserted in the box at position /, 
after which a measurement is performed to find m particles in the left room. As long as the 
number of particles on each side are unequal, the wall would move to equilibrium position 
and perform work through isothermal expansion. Finally, the wall is removed and the whole 
thermodynamic cycle is completed. All the processes of the cycle are isothermal quasistatic 
processes. 

Szilard analyzed this engine and concluded that the entropy produced during the mea¬ 
surement saves the second law of thermodynamics[? ]. Bennett and Landauer argued that 
erasuring information (demon’s memory) generates entropy ks In 2, where ks is the Boltz¬ 
mann constant, which compensates the entropy decreasing in the thermodynamic cycle[? 
? ]. Now it is widely accepted that the total entropy generated during the measurement 
process and the erasure of demon’s memory compensates the entropy decrease of the engine 
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[? ].This reveals the relationship between physics and information, which leads to much 
interest in information in thermodynamic cycle and quantum heat engine [???]. 

Recently, Kim et al. gave a formula of the total work in a thermodynamic cycle for the 
quantum Szilard engine (QSZE), which is related to the relative entropy in the classical limit 
[? ]. In the QSZE with two indistinguishable identical particles in the infinite potential 
well they found that indistinguishability plays an important role in the total work in a 
thermodynamic cycle [? ]. Lu et al. generalized this QSZE to contain arbitrary number of 
identical particles [? ]. They found that in high-temperature limit, the works in a cycle are 
equal for boson and fermion QSZE, which can be understood by the transition from quantum 
to classical statistics when temperature increases. However, in low-temperature limit, the 
physical scenario between fermion and boson QSZE is quite different. The fermion QSZE 
shows the parity effect, namely, the QSZE can extract work in a thermodynamic cycle for 
odd number of fermions, but cannot extract work for even number of fermions. For boson 
QSZEs, there exists a phase transition at a critical temperature for extracted work to be 
positive or negative. [? ] These effects of QSZE imply the significant role of spin nature in 
low-temperature limit. However, QSZE studied by Kim’s and Lu’s group ignore the spin 
degree of freedom of particles. In fact, it can be seen later that the spin would result in 
the degeneracy of energy levels, which would alter the configuration states of the system 
in low-temperature limit. So it should be important and interesting for the spin effect in 
QSZE. 

In this paper, we will generalize QSZE to include arbitrary spin and arbitrary number of 
fermion and boson working substance. We will give the analytic expressions of the extracting 
work in a thermodynamic cycle in low-temperature approximation, in terms of the number, 
spin and temperature. We will also show a detail analysis of the thermodynamic properties 
in a thermodynamic cycle with the particle number, temperature and spin of fermion and 
boson, including the total work, absorbed work and phase diagram of spin QSZE (SQSZE). 
We introduce a new efficiency and give the condition of maximizing work of SQSZE. 


II. MODEL OF SPIN QUANTUM SZILARD ENGINES 

We consider a spin QSZE (SQSZE) that contains N fermions or bosons with spin s as the 
working substance working in a one-dimensional infinite potential well. The thermodynamic 
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cycle of SQSZE consists of four processes, namely, insertion, measurement, expansion and 
removal for four distinct states. The total work in a thermodynamic cycle can be expressed 
as [? ] 

= ( 1 ) 
m =0 V JmJ 

where f m = is the probability for having m particles on the left side after measurement 
and Z(£) = J2m =o ^m(^) represents the partition function for the whole system before mea¬ 
surement. Z m (£) is the partition function for the case in which the wall is at the position i 

7 , (£ rn ) 

and m particles are on the left and N — m on the right. f * n is defined as f* n = ; wh ere 

= J2n=o Zn(£Tq) an d Z n (£™) represents the partition function for the case that the 
wall is at the equilibrium position and n particles are on the left, ks is the Boltzmann 
constant and T is the temperature of the heat bath. 

The partition function can be written as Z — Y2 S e~ Ea ^ kBT where the sum runs all states 
of the system. However, in low-temperature limit, i.e. T —>■ 0, it can be assumed that all 
particles occupy the possible lowest-energy states and all systems in the ensemble are in the 
ground state. Namely, the high energy levels and excited states of the system can be ignored 
in the low-temperature limit. The partition function can be rewritten as 

Z = g G e- EJkBT (2) 

where g G is the degeneracy of ground state of the system and E g is the ground-state energy 
of the system. 

Without losing generality, the wall is assumed to be inserted at the position i — J ~. 
where L is the well width, and initially particles are working in the energy eigenstates, 

2 2 fc2 

E n {£) = ; where n = 1,2...labels eigenenergy levels. 

A. Spin fermion case 

Let us consider N spin-s fermions as the working particles. For convenience, we define 
variable u = s + l/2so that u is an integer. Suppose that the total particle number is 
N = 4im + k, where 0 < k < 4 u. The fermions occupy the energy levels following the 
Pauli’s principle. Namely, each energy level of infinite quantum well can be occupied with 
at most 2s + 1 particles. In general, there are two cases for k, 0 < k < 2 u and 2 u < k < 4 u. 
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FIG. 2. The energy state of the fermion SQSZE when the wall is inserted at position L/2 in the 
low-temperature limit 


They have different configuration states and corresponding different partition functions. We 
will discuss these two cases respectively in the following sections. 

Case 1: 0 < k < 2 u, suppose that m particles are on the left side of the well after 
measurement. We can conclude that in the low-temperature limit it is impossible for 0 < 
m < 2 un or 2 un + k < m < Aun + k , since Z m (L/2) -C Z{L/2) and f m 0. So for the 
case that the wall is at the initial position l = L/ 2, there are 2 un particles occupying the 
energy levels from 1, 2,... to nth states on the left and right sides of the well respectively, 
and k particles occupy the (n + l)th state on the left or right side of the well, as is shown 
in Fig. 2. It implies p — m — 2un particles occupying the (n + l)th energy level of the left 
well where p is within 0 < p < k. By using the combinatorics to calculate the degeneracy 
of the ground state (see Appendix A), and the partition function, we can obtain 


h 


2 W 2 " , 0 <p < k 

0 , otherwise 


( 3 ) 


in the low-temperature limit, where C b a represents combinatorial number. 

When the wall moves to the position Z™, there are 2 un particles occupying the energy 
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levels from 1,2,... to nth states on the left and right sides of the well respectively. For the 
case that m particles are limited on the left side, p particles occupy the (n + l)th state on 
the left side of the well and k — p particles occupy the (n+ l)th state on the right side of the 
well. For the case that particles are not limited on one side dne to the tunnelling effect [? 
], remain k particles occupy the (n + l)th state on one side of the well, in which the energy 
of the (n + l)th energy level is lower than that of the other side. So for k is odd, we get 

(-1V syk—p 

f m = 2 " (4a) 

^2 U 


and for k is even, 


/ * 
m 


r<p r ik ~p 

'-'0.11 '-'Om 



e -HpA E ™ +1 



(4b) 


where A E™ +1 = | E n+ x — E n+1 (L — |, p — min {p, k — p} denotes the smaller of the 

number of particles on the (n + l)th level of both sides. 

In particular, when the wall reaches the end of the well, n = 0 and p = 0 or k, f * n = 1. 
By combining the Eqs.(l)(3) and (4) we can obtain the total work in a thermodynamic 
cycle of SQSZE, 


W tot 


k B T In ( g|- 


-2E2p^AO 


for k = odd 


k B T 


1 - 


[c^y 

/^k 

u 


In 


for k 


( 5 ) 


= even 


in which we agree that the sum equals to zero if its upper limit is smaller than its lower 
limit and we obey this agreement in this paper. 

Case 2:2 u < k < 4u, in this case, notice that there are more than 2 u particles occupying 
the (n + l)th energy level in the left or right well, which means that particles are more than 
unoccupied states, we use the hole representation of the state to calculate the partition 
function for convenience. There are q = 2 u(n + 1) — m unoccupied states on the (n + l)th 
energy level of the left well, where q is within 0 < q < 4u — k. By the similar combinatorics, 
(see Appendix A), we obtain 


frr 


, 0 < q < 4u — k 
, otherwise 
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( 6 ) 



In the same way, we get 


for odd k and 
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4 u — k 
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c: 
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, q^ 
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4u—k 

2 

4 u—k 
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(7a) 


(7b) 


for even k, where A E™ +1 = \E n+1 — E n+1 (L — i™') |, A = min { 4 m — q — k, q} denotes 

the smaller of the number of unoccupied states on the (n + l)th energy level of both sides. 
By the same way, the total work in a thermodynamic cycle of SQSZE can be expressed 


as 


W tot = 


kuT In 
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for k = odd 
for k = even 


( 8 ) 

Notice that the total works have the same form for both cases, we can rewrite them as 

Wt 0 t = DpksT — Wqp ( 9 ) 

where Dp is proportional to the change rate of total work as temperature increases and Wqf 
represents the absorbed work of fermion SQSZE at zero temperature. 

For 0 < k < 2 u, 


Dp = 


1 

kp 



and for 2u < k < 4m, 
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for k = odd 
for k = even 


D F = 


1 

kp 



In 


c; 


c: 


4u 

4 u — k 


1 - 


c. 


(4u-k)/2\ 2 


ct 


In 


cl 


cl 


for k = odd 
for k = even 


( 10 ) 


( 11 ) 


B. Spin boson case 

For boson SQSZE, suppose that N spin-s bosons are working in SQSZE and m of them 
are in the left well after measurement. Each particle has 2s + 1 different spin states, which 
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leads to 2s +1 degenerate degree of freedom on each energy level. By the combinatorics (see 
appendix B), in the low temperature approximation, we can obtain 


fin 


s (~i2 s 

^ m.+2s'- / N-m.+2s 
/^i4s+l 
°Af+4s+l 


( 12 ) 


Similarly, we get 


/ '* 
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1 


m/0 and m ^ A 
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(13a) 


for odd A, and 
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N 
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(13b) 


for even A, where A E™ = \E\ (7™) — E\(L — lf q )\, /J = min {m,N — m} denotes the 
smaller of the number of particles on the first level of both sides. 

Thus we can obtain the total work in a thermodynamic cycle in SQSZE 
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(14) 


Similar to the fermion case, the total work in Eq. (14) can also be written as 


W tot — DsksT — Wqb 


(15) 


where 
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(16) 


is proportional to the change rate of total work as temperature increases and Wqb represents 
the absorbed work of boson SQSZE in the low temperature approximation. 




FIG. 3. (Color online) The change of Dp and Wq F along with N, in which u= 5, L = 10 9 m, 
M = 10~ 26 kg 

III. BASIC PROPERTIES OF D F ,W 0F ,D B AND W 0B 

A. Fermion SQSZE 

For the given width of the quantum well L = lnm, the mass of particle M = 10~ 26 kg, and 
the particle spin s = 9/2, we investigate D F and Wq F versus N, which is shown in Fig.3. It 
can be seen that D F , which represents the change rate of total work as temperature increases, 
shows even-odd oscillation and periodicity, which modifies the parity effect of the spinless 
fermion QSZE[? ]. Besides, there also exists absorbed work in fermion SQSZE, which does 
not happen in spinless fermion QSZE. These behaviors of Wq F and D F imply quantum 
nature of fermion SQSZE in low temperature. In high temperature, both periodicity and 
even-odd oscillation of Wq F and D F will vanish)? ]. 

Notice that the absorbed work increases as the particle number increases, we define the 
average absorbed work by 




W 0F 

N 


( 17 ) 


We investigate the average absorbed work versus the number of particles in the Fig.4(a). 
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FIG. 4. (Color online) (a) The change of Wqf along with N. (b) The partial enlarged diagram of 
(a) when N is very big. u= 5, L = 10 _9 m, M = 10 26 

The Wqf declines rapidly with oscillation in the range of small number particles. As the 
number of particles increases, the Wqf becomes periodic oscillation and trends to exact 
periodicity in the large-number limit of particles, as shown in Fig.4(b). (See Appendix 
C) This implies that the average absorbed work of fermion SQSZE is neither an intensive 
quantity nor an extensive quantity, which is also an unique characteristics of spin ferrni 
system. 

We can also investigate the maximum of D F and Wqf versus spin u = s + l/2of particles 
shown in Fig. 5. The maximum of Dp and Wqf increases as spin of particles grows. It 
implies that the spin of particles plays an important role in the total work of fermion SQSZE, 
which is also a quantum nature of of fermion SQSZE. 

As shown in Fig.6, we can explicitly plot the relationships between the total work W tot , 
temperature T and the particle number N. It can be seen that the total works of the 
different-spin fermion SQSZE oscillate with the particle number increasing, and increase 
with temperature. For s = 1/2 case, the total work W to t > 0 in the whole T — N region, 
which implies that there is no absorbed work. For s > 1/2 cases, the total works can be 
negative in the large-N regions, which can be viewed as a phase transition induced by the 
absorbed work variation in the high-spin fermion SQSZE, which can be also seen in the 
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FIG. 5. (Color online) (a) The change of the maximum of Dp along with u(u = s + 1/2). (b) 
The change of the maximum of Wqf along with u when n has different values. u= 5, L = 10 -9 m, 
M = 10 ~ 26 kg 

phase diagram in Fig. 9. 

B. Boson SQSZE 

For boson SQSZE, the Dp and Wqb versus N for given different spin s is shown in Fig. 7. 
It can be seen that both Dp and Wqp increase with the particle number, and there is some 
oscillations when the spin increases. Interestingly, Dp increase with spin increasing, but 
Wqp decreases with spin increasing. These properties reveal the difference between boson 
and fermion SQSZE. 

In the large spin limit, Dp and Wop can be expressed as 

| IV In 2 for N is odd 

r N/2\ (18) 

1- —pp- i N In 2 for N is even 

EJTi mC N' f o r N is odd 

hm W 0B = < iv_ (19) 

s ^°° y Em= i rnC^AE] 11 for N is even 

To explicitly understand the relationships between the total work W to t , temperature T 

and particle number N, we plot the Fig.8. It can be found that the total work W tot of the 
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FIG. 6. (Color online) The change of total work Wtot along with temperature T and particle 
number N when (a) u= 1, (b)u=2, (c)u=5, (d)it=10. L = 10 ~ 9 m, M = 10 ~ 25 kg 



FIG. 7. (Color online) The change of Db and Wqb along with N when s gets different values, in 
which L = 10~ 9 m, M = 10 ~ 26 kg 

boson SQSZE increases with temperature and spin increasing, but has no great oscillation 
with the particle number. There is a low-temperature region, in which the total work is 
negative. As spin increases, this region becomes small, which can be seen clearly from the 
phase diagram in Fig. 9. The basic behavior of the total work of the boson SQSZE is quite 
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FIG. 8. (Color online) The change of total work Wtot along with temperature T and particle 
number N when s=0,l,2,5. L = 10 ~ 9 m, M = 10 ~ 2e kg 

different from fcrmi SQSZE, especially for high-spin working matters. 

C. Phase diagram of SQSZE 

As a heat engine, an important issue is the relationship between work, temperature and 
particle number. We define the critical temperature of fermion SQSZE to study the working 
properties of SQSZE, 


Wof 

T cF = where N + Aun (20) 

Dpk,B 

We numerically investigate the working properties of fermion SQSZE in terms of the 
critical temperature and the particle number in Fig. 9, which can be regarded as a phase 
diagram of fermion SQSZE. The curve in the phase diagram goes up and shows a periodic 
oscillation. The upper region of the curve in the phase diagram denotes the fermion SQSZE 
engine doing positive work, while the lower region of the curve represents the fermion SQSZE 
doing negative work. 

Similarly, the critical temperature of the boson SQSZE can also be defined by 

T cB — -yp-y—, where N ^ 0 (21) 

Db^b 
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FIG. 9. (Color online) The work phase diagram of the fermion QSZE, in which u= 5, L = lCU 9 m, 
M = 10 _26 %. 

The phase diagram of the boson SQSZE is plotted in Fig. 10 for given spin of particles. 
The critical curves for IT = 0 go up with oscillation and separate for different spins as 
the number of particles increases, but it looks quite different from the curve of the fermion 
SQSZE. 

Actually, Lu and Long found that the phase transition occurs in the spinlcss boson 
SQSZE[? ], and there have been some works discussing the superposition effect in a two- 
state quantum engine [? ? ], its effeciency[? ? ]. However, all of these works ignore the spin 
effect in the quantum engine. Here, we found that the spin degree of freedom can decrease 
the critical temperature of the phase transition in boson SQSZE. This is because the spin 
degree of freedom will modify the configuration states of the working matter. 


IV. INFORMATION-ERASURE WORK AND INFORMATION-WORK EFFI¬ 
CIENCY 

According to the Landauer’s erasure principle, the erasure or resetting of memory re¬ 
sults in the energy dissipation, which would cause entropy increase. The erasure work is 
proportional to the amount of information stored in memory, [? ] 
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FIG. 10. (Color online) The work phase diagram of the boson QSZE when s gets different values, 
in which L = 10 ~ 9 m, M = I0~ 2e kg. 


W eras = k B T\n2H(P) (22) 

where H(P ) = — X^=i -f* log 2 Pi represents Shannon entropy (or information entropy) in 
unit Bit. Due to the quantum indistinguishability of particles, the demon does not know 
where particles are. In fact, the information that demon needs to acquire in the measurement 
process is the number of identical particles on each side. Since f m is the probability of finding 
m particles on the left side, the erasure work is given by 

N N 

Weras = ~k B T In 2 Y f m log 2 f m = ~k B T ^ f m In f m (23) 

m =0 m =0 

Suppose that the erasure and resetting of the memory happen in every thermodynamic 
cycle, and the erasure work is provided by the engine. The net work in a thermodynamic 
cycle of SQSZE can do is given by 

N 

w net = W tot - W erCLS = k B T Y, fm In r m (24) 

m=0 

It can be seen that the profound relationship between the net work in Eq.(24) and the 
second law of thermodynamics. Since ff < 1, the net work done by SQSZE must be 
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negative or zero, which corresponds with the Kelvin’s statement of the second law of ther¬ 
modynamics because the heat engine cannot extract any work from a single heat source in a 
thermodynamic cycle, after resetting the whole system including the memory to the initial 
state. This irreversibility can also be understood by the tunnelling effect, which leads to the 
irreversibility as diffusion does. When f* n — 1 for given m , i.e Z m (l™) = Z(l™), there is no 
tunnelling or tunneling does not take effect for the case that m particles are on the left side, 
so unless all the f* n equal to 1, the SQSZE must be irreversible due to the tunnelling effect. 

We hope that the work done by SQSZE is as much as possible while the erasure work 
is as little as possible. In order to judge the quality of a SQSZE from this aspect, we 
define an efficiency as the ratio of the total work and the erasure work, we called it as the 
information-work efficiency. 


_ W tot 

r,1 ~ w w 

v y erat 

Notice that the total work done by QSZE might be negative, the domain of information- 
work efficiency is (—oo, 1]. When the information-work efficiency is negative, it means that 
the engine cannot do any work. When the information-work efficiency is positive, it implies 
the engine can do work even though it still needs erasure work. When the information-work 
efficiency is higher, we believe that the engine is more efficient and more reversible because 
less energy is dissipated. In particular, the engine become reversible when the information- 
work efficiency reaches 1 because there is no energy dissipation. 

Substitute Eqs.(3)(6) (9) and(23) into Eq.(25), we obtain the information-work efficiency 
in fermion system 



Vf-v 
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_ hT’Ep=o/p' n /p 

DphsT—WoF 


for 0 < k < 2u 
for 2u < k < Au 


(26) 


(-rP (~<k-V 

where f p = ~ 2 % 2u and f q 
system is given by 
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Likewise, the information-work efficiency in boson 
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V. THE CONDITION OF MAXIMUM EXTRACTED WORK AND HIGHEST 


INFORMATION-WORK EFFICIENCY 


Naturally, we want to know the condition of maximum extracted work and highest 
information-work efficiency. From Eqs.(9) and (15) we know that the extracted work can 
reach maximum when the engine does not absorb work from environment, namely Wqf = 0 
(or Wob = 0). It can be verified that Wqf = 0 only when k — 0,1,2,4 u — 2,4m — 1 and 
Wqb = 0 when N = 0,1, 2. So we can verify that, for k — 1, 4 m — 1 the maximum extracted 
work of fermion SQSZE is W^f 10 " = k B T In 2, and for N = 2 the maximum extracted work 

B,max _ r-p 2s+2 | n 4s+3 


of boson SQSZE is W t n J nax = k B T In Interestingly, W t 


B,max 

tot 


= AmT§4§ In > 


4s+3 s+1 

k B T In 2, which can be generalized to W t ^ max > W t ^( na: ' for all cases. This means that we 
should use bosons as working substances in order to obtain maximum extracted work. 

For the information-work efficiency defined by Eq.(25), it is possible to reach maximum 
when the engine does positive work, i.e. k = 0,1, 2, 4 m — 2, 4 m — 1 or N = 0,1, 2. By using the 
Eqs. (26), for fermion SQSZE, the maximum information-work efficiency can be obtained 
vf-™* — 1 for k — 1,4m — 1. Similarly, for boson SQSZE, the information-work efficiency 
reaches maximum, r]fj_™ ax = 1 for N = 1. 

Interestingly, the second-highest information-work efficiency of both fermion and boson 
SQSZEs can be expressed as a same form, 

1 


F(B),h2 = 
li—w 


1 ! (1—«f(s)) ln(l-a_F( S )) 

1 ' i “F(B) 

a F(B ) In Hi 


(28) 


where = ^-j- with k = 2 or 4 m — 2 for fermion SQSZE, and a B = fyp| with 

N = 2 for boson SQSZE. 

Therefore, the working particle number plays an important role in the extracted work 
and information-work efficiency for SQSZE. If we can control the particle number in 4nn +1 
or 4nn + 4 m — 1, where n e N, the fermion SQSZE has the maximum extracted work 
and highest information-work efficiency. However, for the boson SQSZE, when the working 
particle number is N = 2, the boson SQSZE has maximum extracted work and the highest 
information-work efficiency corresponds to the working particle number N = 1. Moreover, 

when rji _ w = 1, the thermodynamic cycle is reversible, namely, when N = Aim + 1, or 

4 un + 4m — 1 the fermi SQSZE and N = 1 for the boson SQSZE, the QSZE becomes 
reversible. 
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VI. CONCLUSION 


The QSZE plays an important role in understanding the fundamental physics of the 
quantum engine, information, and quantum thermodynamics. We generalize QSZE to spin 
QSZE (SQSZE) and systematically study the basic physical properties of both fermion and 
boson SQSZEs in low temperature approximation. We give the analytic formulation of the 
total work, information erasure work, and efficiency. Based on this formulation, we find a 
series of the physical properties of SQSZE: 

(1) The total work depends on both particle number and particle spin. The change rate 
Dp as temperature increases exhibits periodicity and even-odd oscillation with the particle 
number. The parity effect embeds in the periodicity and even-odd oscillation. This is a 
generalization of spinless fermion QSZE, in which the total work depends on only particle 
parity. [? ] 

(2) For the fermion SQSZE, the average absorbed work from environment osicillates regularly 
and periodically in large-number limit, which implies that the average absorbed work is 
neither an intensive quantity nor an extensive quantity. This is a novel and pure quantum 
phenomenon in quantum thermodynamics. 

(3) The phase diagrams of both fermion and boson SQSZEs gives the SQSZE doing positive 
or negative work in the parameter space of the critical temperature and the particle number 
of system. The critical temperature of phase transition depends sensitively on the particle 
number. 

(4) By using the Landauer’s erasure principle, we give the erasure work in a thermodynamic 
cycle and define a new efficiency (we refer it information-work efficiency) to measure the 
engine’s ability of utilizing information to extract work. 

(5) We give the particle number conditions under which the maximum extracted work and 
highest information-work efficiencies for fermion and boson SQSZEs can be achieved. 


These results reveal the basic physics of SQSZE, especially for characteristic difference 
between spinless QSZE and spin QSZE, as well as fermion and boson working particles in 
SQSZE. These physical properties of SQSZE provide some hints to understand the relation¬ 
ships between quantum mechanics, quantum thermodynamics and information physics. 
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Appendix A: DERIVATIONS OF f m and f* m IN FERMI SYSTEM 


Let us consider N fermions with spin quantum number s occupying the energy levels. 
Defining u = s + | and suppose that N = Aim + k where n > 0, we know that each spin-s 
fermion has 2u different spin states, which is equivalent to 2 u degrees of degeneracy on each 
energy level. 

Derivation of f m 

In low-temperature approximation, we ignore the temperature fluctuation and consider 
only the ground state in the canonical ensemble, namely, particles occupy as low energy 
levels as possible and are constrained by Pauli’s exclusive principle. Suppose that the engine 
system is an infinite potential well and in the first step of SQSZE, a wall is inserted at a 
half width position. Thus, there are 4 un particles occupying the low energy levels 1,2, ...n 
in the left and right sides of well, and the other k particles occupy the (n + l)th energy level 
in the left or right side. Namely, k particles occupy two energy levels with 2s + 1 degrees of 
degeneracy . Using the combinatorics, the number of configuration states is C± u . Therefore, 
the partition function Z(- 1) in the low-temperature approximation is given by 



where E a (|f) = is the partition function for m particles on the left side. 

For 0 < k < 2 u, the domain of m is 2 un < m < 2 un + k, and for 2 u < k < 4 u, the domain 
of m is 2im + k — 2u<m< 2un + 2u, since if m equals to other values, Z m (L/2) <C Z(L/2) 
and f m —>• 0, which means it is impossible for m particles located on the left side. Suppose 
that p — m — 2un particles occupy the (n + l)th energy level on the left side, namely k — p 
particles occupy the (n + l)th energy level on the right side of the well, the partition function 
Z m ( f) is given by 


Z 


m 




(A2) 


Notice that for 2 u < k < 4 u, it is convenient to express the number of configuration states 
in terms of the hole-representation, we obtain the probability 
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( /~iv r ik ~P 


fm = 


< 2 ) 


for 0 < k < 2u and 0 < p < k 


Z(f) 


= < 


°4 u 
,~<q f-iAu — k — q 

2u Ju~k — f° r 2u < k < Au and 0 < q < Au — k 


c: 


(A3) 


4 u 


0 for otherwise 

where q = 2u — p = 2u(n + 1) — m represents the number of unoccupied states on the 
{n + l)th energy level of the left side. 

Derivation of /*, 

For f* z , the wall expands to a balance position £™. It satisfies the balance equation [? ? 


E F »(c a ) 


n= 1 


dEn (q) 

SIT, 


E ( L - O 


dE n (L - C a ) 


(A4) 


n =1 


d(L- e%) 

where P n (£™ q ) represents the average number of particles on the nth energy level on the left 
side, and P n (L — £™ q ) on the right side. In the low-temperature approximation, particles 
occupy as low energy levels as possible. Thus, the Eq. (A4) can be reduced to 


fm 

L eq 


2«Efch 2 + ([£] + i) 2 « 


(A5) 

2« EU- '« 2 +([^]+l) 3 6 

where [x] = max(i e Z, i < x), and a = mod (m, 2 u) and b = mod (N — m, 2 u). The Eq. 
(A5) can be simplified to 


fm 
^eq 


un(2n + 1) + 3 p(n + 1) 


L — £™J un(2n + 1) + 3 (A: — p)(n + 1) ^ 

It can be verified that E n+1 (£™) > E n (L — £™), and E n (£™) < E n+1 (L — £™), so the 
2 n lowest degenerate energy levels of the systems are Ei(£™)E 2 (£™) ■ ■ ■ E n (£™) and E\(L — 
£™)E 2 (L-£™)---E n (L-£™). 

Case (1): 0 < k < 2u 

When k is odd and p < |, the partition function of the system for the wall at the balance 


-eq 


IS 


Z = ^ e -/3[ 2 «( Ei=i E Ar q )+E7=i E iP-^ q ))+ kE ^+i( L -^ q )] 
and similarly, the partition function for m particles on the left side is obtained 

Zm = (jP^(jk-P e -P[2u(j2?=i E i(t?q)+E?=iEi(L-e™ q ))+pEn + i{e™ q )+(k-p)E„ +1 (L-e™)] 


(AT) 


(AS) 


20 



Consequently, the probability is obtained 


f* = Zm (^) = C 2u C 2u P -p0(E n+1 (tZ,)-E n+ i(L-t? n )) 

Im Z (£*j) C* 


(A9) 


When k is even and p > |, similarly to the Eq. (A9), we obtain 


/ '* 
m 


2u L '2u -( k-y)t3(E n+1 (L-i™)-E n+1 (l”)) 

r<k c 
°2 u 


Combining the Eq. (A9) and Eq. (A10), the can be rewritten as 


(A10) 


r<p r k ~p 

f* = -pftA g”, 

J m k c 

°2u 


(All) 


where AE™ +1 = |A n+1 — E n+1 (L — |, p = min {p, k — p} denotes the smaller one 

of the number of particles on the (n + l)th level of the both sides. 

However, for 0 < k < 2 u, and k is even, p = |, = L — = A the (n + l)th energy 

level of both sides is the same, we get 


Z (l m ) — C k p~ /3 [ 4u ^i=l Ei (^) +kEn+1 (^)\ 

V eq) 4 u 

(A12) 

and 


Z m (l™) = c\ u c\ u e-^ u ^= i^(#)+fc^+i(#)] 

(A 13) 

Consequently, we obtain 

k k 


r<2 r<2 

r* _ U 2w L/ 2m 

Jm ~ s-ik 

°4 u 

(A 14) 

Case (2): 2u < k < Au 

Similar to the case (1), when k is odd, we get 


/~<q ^iu—q—k 

f* _ L/ 2n L/ 2« -A/3AEr + i 

Jm s~i4u—k 

°2 u 

(A 15) 

and when k is even and q ^ 4u ~ k , we have 


/~iq s~iiu—q—k 

f* _ L/ 2u L/ 2u -A/3AEr + i 

Jm s~iAu—k 

°2 u 

(A16a) 

When k is even, but q = 


4u — fc 4u — k 


2 2 

r* _ Ky 2u ^2u 

Jm s~iAu—k 

°4 u 

(A16b) 
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where A = min {4 u — q — k,q} denotes the smaller of the number of unoccupied states on 
the (n + l)th level of the both sides. 

In particular, when n — 0 and p — 0, k, the wall is moved to the end of box so f* n = 1. 


Appendix B: DERIVATIONS OF f m and fa IN BOSE SYSTEM 


Suppose that the system has N bosons with spin quantum number s, we know that 
each boson has 2s + 1 different spin states, which leads to 2s + 1 degrees of degeneracy on 
each energy level. In the low temperature approximation, we also ignore the temperature 
fluctuation and consider only the ground state in the canonical ensemble. 

Derivation of f m 

Notice that when the number of particles is a and the degrees of degeneracy is w, the 
number of micro-states is given by = C“ +a _ 1 [? ], we get 


Z 



/'~y4s+1 

'~ / V+4s+l 


e 


4V/?Ei(f) 


(Bl) 


and 



f<2s -NpE\(- 1 ) 

( - / m+2s'- / jV-m+2s e 


(B2) 


Therefore, we obtain 


fm. = 


z m (f) 

z(f) 


s r<2s 

'- y m+2s'- y JV— m+2s 
<^y4s+l 
°Af+4s+l 


(B3) 


Derivation of f* n 

In the boson case, the Eq. (A4) becomes 


fim 
^eq 


trm 
^eq 


m 


N 


m 


(B4) 


Since the lowest energy level of boson system is Ei(£™) or E\ (L — £™). When N is odd 
and 0 < m < we get 


z K) 


r N NpE^L-P™) 

°Af+2s e V ’ 


(B5) 


and 


(q) 


q 2 , -0[mE 1 (<3)+(W-m)E,(l-<3)] 

777.—|— Jff 7T2-—|— 2 S 


(B6) 
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Consequently, we obtain 


te) 

z («?,) 


^'m+2s^'N—m+2s — mf)\E\(i ™ q ) — E\ (L—l™)] 


riN 
Ly N+2s 


(B7) 


Similarly,when IV is odd and A±I < m < N, we can obtain 


f* _ ^m+2s^-m+2s (N-m)/3[E 1 (L-e™)-E 1 (e™)] 

Jm /~iN e 

°7V+2s 

Combining the Eqs. (B7) and (B8), we get 


(B8) 


where A E™ 


/ * 
m 


r^2 s /^i2s 

^ m+2s^ N-m+2s -iiffAE? 1 

c*n e 

Ly N+2s 


(B9) 


Ei (£™) — Ei (L — , /r = min {m, N — m} denotes the smaller of the 


number of particles on the Erst level of the both sides. 


By the same way, when N is even and m ^ y, 0, A, we get 


/~i2s (~i2 s 

r* ( - y m+2s'~'N—m+2s — fj,pAE™ 

Jm ~ N e 

Ly N+2s 


(Bio) 


However, when N is even and m — y, l™ q — L — — |f, the Erst energy level of both 


sides is the same, so we get 


In particular, when m 


/ '* 
m 


c 


2s /~<2s 
f+2s L 'f+2s 

r<N 

C/ V+4s+l 


(BH) 


0 or m — N, the wall reaches the end of the box, which makes 


f* = 1 . 

J m 


Appendix C: DERIVATION OF lim^^oo Wqf 


For the infinite potential well model, A E™ + i is written as 


AK+i = 


(n + l) 2 7 r 2 h 2 


2 M 


te ) 2 (l - c ,) 2 


It can be rewritten as 

(n + l) 2 Eh 2 


AK+1 = 


1 + 


L-IV 


2 ML 2 


L-i? q 


1 + 


r _/m 
^ L eq 


+ 




jm 
**eq 


L-l™ 


(Cl) 


(C2) 
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By substituting Eq. (A6) to Eq. (C2), we get 


A K 


n +1 


4tt 2 H 2 
3 ML 2 


(n 


l) 2 


3 (k - 2 p) (n + 1) 


tin (2 n + 1) + 3 (k — p) (■ n + 1) 


When n> 1, the equation (C3) can be approximated to 


(C3) 


A-ji 2 K 2 


3 ML 2 
Case (1): 0 < k < 2u\ 


3 (k - 2i) 

4vr 2 K 2 

( n + 1) 

k p 

un(2n+l) , 3 (k—p) 

(n+1) 2 + (n+1) 

" ML 2 

2 u u 


When k is odd, in the large n limit, W 0 p is given by 

(A _ E) 2fc : 

V 2 u uJ 7T fl 


v 7 TT~ 87r2 ^ 2 V 

lim Wqf = ——— lim 

n —>oo JVl Jj* n—>oo 


_ /~<p 

2 ^ °2u°2u 

p=0 P 


r~tk 

4im+fc 

n+1 


fc-i 

2fc2 2 x-yp 

'- y 2u'~ y 2u 


(C4) 


u 2 M L 2 


5>- 

p=0 


r<k 

°4n 


(A; — 2p) (C5a) 


Similarly, when k is even, we have 


lim Wqf 

n—>• oo 

Case (2): 2u <k < Au 


71-2 ft 2 \ ^ C V 2yp2u 

u 2 ML 2 ^ Cl 


(.k — 2 p) 


In the same way, when k is odd, we have 


(C5b) 


lim II'op 

n—> oo 


7 T 2 K 2 
u 2 ML 2 


4 u — k — 1 
2 

E * 

9=0 


piq f~i4u—k—q 
Ly 2u Ly 2u 
s~i4u—k 
°4 u 


(4m 


k — 2 q) 


When k is even, we get 


(C5c) 


IT 


■h 2 


--i 


}™ W °r= u 2MF ^ 

q =o 


4 u — k 
2 

\ " Ky 2u KJ 2u 


f~iq ^4u—k—q 
Oo,, Oq 


Qn 


4n—fc 


(4m — 


We can prove that ITq/t becomes periodic and symmetric about k 


k - 2q) (C5d) 

= 2 u in the large n limit. 
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